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The universality class, even the order of the transition, of the two-dimensional Ising model depends 
on the range and the symmetry of the interactions (Onsager model, Baxter- Wu model, Turban 
model, etc.), but the critical temperature is generally the same due to self-duality. Here we consider 
a sudden change in the form of the interaction and study the nonequilibrium critical dynamical 
properties of the nearest-neighbor model. The relaxation of the magnetization and the decay of the 
autocorrelation function are found to display a power law behavior with characteristic exponents 
that depend on the universality class of the initial state. 



I. INTRODUCTION 

The parameters of a physical system are often subject 
to a sudden change, as for example fast cooling or heating 
or the switching on or off of an external field. In material 
science these processes are used to prepare new types of 
(glassy) states of matter. After such a quench the sys- 
tems are out of equilibrium and their dynamical proper- 
ties can be completely different from those known in the 
equilibrium situation. The phenomena of aging and reju- 
venation are typical examples of such glassy dynamics^. 

In a theoretical approach a simple and often stud- 
ied example 3 -! 4 -! 5 - is a d-dimensional ferromagnet, which 
is quenched from its high-temperature (paramagnetic) 
phase to a temperature T which is below or at the criti- 
cal temperature T c . For T < T c we have phase-ordering 
kinetics where order grows through domain coarsening 6 . 
For a critical quench, T — T c , the domains are fractals, 
and the growth and dynamics involve critical exponents 7 
such as the magnetization scaling dimension x = (3/v (/? 
and v are the critical exponents of the magnetization and 
the correlation length, respectively) and the dynamical 
exponent z. 

In this case one often measures the autocorrelation 
function^, G(t,s) = (a(t)a(s)), where a(t) is the oper- 
ator of the magnetization and t and s are the obser- 
vation time and the waiting time, respectively. Aging 
is manifested by the fact that G{t, s) is nonstationary, 
but instead has the scaling form G(t,s) — t~ 2x / z G(s/t), 
where the scaling function for small argument behaves 
as lim T ^o G(t) ~ T ( d ~ x i- X )/ Z , Here Xi is a new nonequi- 
librium exponent, the anomalous dimension of the initial 
magnetization's. The autocorrelation function for t ^> s 
behaves as Git) ~ t~ x l z , where the autocorrelation expo- 
nent satisfies the relation A = d — Xi + x. Another com- 
mon measurement after a critical quench concerns the 
relaxation of the magnetization 9 , m(t) = (cr(t)), where 
the system is initially prepared in a state with a small 
initial magnetization m^. Here we have asymptotically 
m(t) ~ t , where the initial slip exponent is given by 
8 = (xi — x)/ z. 



In critical nonequilibrium dynamics the initial state 
might contain some kind of correlations. In the case 
of a perfectly correlated initial state the relaxation pro- 
cess involves only equilibrium exponents. Here we have 
m(t) ~ t~ x/z , and, similarly, for t » s, G{t) ~ t~ x/z . 
The dynamic crossover between the ordered and the dis- 
ordered cases has been the subject of a recent series 
of paper s 11 ' 12 ' 13 . Another possibility is given by initial 
states that display quasi-long-range order, i.e. where 
correlations decay as a power law. This happens for 
the two-dimensional XY model if both the initial tem- 
perature, Ti, and the final temperature of the quench, 
T, are below the Kosterlitz-Thouless temperature, Tkt- 
If Ti < T <C Tkt, we have according to spin- wave 
theory^ G(t) ~ t -(*(T)-x(.T z ))/z ^ where x ^ is the 

value of the anomalous dimension at the given tem- 
perature. Other recent studes have focused on the d- 
dimensional spherical model with an initial state of pre- 
scribed correlations^ or on the Ising model with initial 
states generated through random field effects 17 . 

In the present paper we consider quenches during 
which the temperature of the system remains the same, 
but where the form and the local symmetry of the in- 
teractions are changed. Interestingly, the recent progress 
in the experiments on phase transitions in optical lat- 
tices could make this type of investigations possible. We 
here consider the case of Ising spins on a square lattice 
with different types of interactions. The phase transition 
encounterd in the Ising model is considered to be the 
paradigm of an order-disorder transition since the exact 
solution of Onsager— of the standard model with near- 
est neighbor couplings. The same model, however, with 
three-spin product interaction for each elementary trian- 
gle belongs to a different universality class. According to 
the exact treatment by Baxter and Wu 1 ^ the (static) crit- 
ical exponents of this so called Baxter- Wu (BW) model 
are the same as those of the four-state Potts model 20 . 
Still another self-dual model has been introduced by 
Turban^ii^i 2 ^ and other s 24 ' 25 that has nearest-neighbor 
interactions in the vertical direction, but n-spin prod- 
uct interactions in the horizontal direction. Of course, 
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for n = 2 we recover the Onsager problem, whereas for 
n = 3, according to symmetry arguments, approximate 
mappings and numerical investigations, the system be- 
longs to the four-state Potts universality class. For n > 4 
the phase transition is of first order. 

In our study we prepare the system in an initial state 
that is an equilibrium critical state of the BW-model or 
of the multispin model with n = 3 and n — 4. In the ini- 
tial state of the BW-model and that of the n = 3 model 
there are the same type of critical correlations, since the 
two models belong to the same universality class. On 
the other hand for the n = 4 model we have phase co- 
existence at the phase transition point between the or- 
dered and the disordered phases. Having prepared the 
system in this way, we change at time t = the form of 
the interaction and then let the system evolve according 
to spin-flip dynamics with nearest neighbor interactions. 
We thereby measure the relaxation of the magnetization 
and the decay of the autocorrelation function and deter- 
mine the exponents X/z and 9. 

The structure of the paper is the following. In Sec. |TT] 
we define the models and describe their critical proper- 
ties. The results of noncquilibrium relaxation are pre- 
sented in Sees. IIIII and IIVI where the initial state corre- 
sponds to an equilibrium state of a second order or of a 
first order transition point, respectively. We discuss our 
results in the final section. 



II. MODELS AND THEIR CRITICAL 
PROPERTIES 



Hamiltonian is given by 

Ti-BW = — ^ JbW (cijCi-ljCiJ-l + ■ 

The ordered phase of the system is four-fold degenerate, 
with the majority spin orientations in the three equiva- 
lent sublattices being given by : |, |, f; t, i, U U t> li an d 
J., J,, t- According to exact results the critical point of the 
system is located afc^ 

smh(2J BW /kT c ) = 1 , (4) 

and the static critical exponents are the same as for the 
four-state Potts models, but without logarithmic cor- 
rections to scaling. Numerical results indicate that uni- 
versality also holds for the dynamical exponent, z. The 
nonequilibrium exponents, A and 9, however, seem to be 
different, see Table |U 

c. Turban model. In this model we have nearest 
neighbor interactions with strength Ji in the vertical di- 
rection and n-spin product interactions with strength J„ 
in the horizontal direction, so that the Hamiltonian is 
given by 

/ n-l \ 

7~in = - ^ I J^VijCij+X + Jn ]^[ Ci+fe,j I ■ (5) 
i,j \ k=0 I 

The ordered phase of the system is 2™ _1 -fold degenerate. 
The model is self-dua l 21 ' 23 ' 24 and the self-dual point is 
located at 



We consider in the following Ising spins, <x;.j = ±1, 
at the sites of a square lattice with different types of 
ferromagnetic interactions. 

a. Ising model. The standard Ising (or Ising-Lenz) 
model only contains nearest neighbor ferromagnetic cou- 
plings, so that the Hamiltonian is given by 

Ki = ~ / ] J (o-i,j<Ti,j+i + a hj a l+1J ) . (1) 

where i and j label the lattice sites, whereas J is the 
strength of the couplings. The critical point is given by 
the condition^ 

sinh(2J/fcT c ) = 1 , (2) 

which separates a two-fold degenerate ordered phase from 
a paramagnetic phase. The static critical exponents are 
known exactly, whereas the dynamical exponent z and 
the nonequilibrium exponents X/z and 9 are calculated 
numerically with high precision. The values of these ex- 
ponents are collected in Table U 

b. Baxter- Wu model. In the Baxter- Wu model we 
have three-spin product interactions with strength Jbw 
between spins located on elementary triangles, and the 



sinh(2 J 2 /kT c ) sinh(2 J n /kT c ) = 1 . (6) 

It is known from numerical studies that a single phase 
transition takes place in the system, thus the phase tran- 
sition temperature coincides with the self-duality point. 
In the following we take J2 = J n {= J = Jbw), so that 
all models described in this Section have the same crit- 
ical temperature, see Eqs.(0), Q and ©. As already 
mentioned in the introduction the Turban model with 
n = 3 has a continuous phase transition which belongs 
to the (static) universality class of the four state Potts 
mode l 22 i 25 ' 26 i 27 ' 28 i 29 , even the logarithmic corrections are 
expected to be of the same form for the two models. The 
numerical estimates of the dynamical exponent of the two 
models are somewhat different, although they could be 
the same within the error of the calculation. The same 
conclusion holds also for the noncquilibrium exponents, 
A and 9, see Table HJ 

The model for n > 4 has a first order transitio n 30 ' 31 . 
Detailed numerical studies are available for the n = 4 
case which exhibits a latent heat of A/kT c = 0.146(3) 
and a jump of the magnetization from zero to m c = 
0.769(6). Nonequilibrium relaxation studies of this model 
have been performed recently^ 2 -. The autocorrelation 
function has thereby been found to approach its limiting 
value, given by the magnetization in the ordered phase 
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TABLE I: Upper part: static and dynamic critical quanti- 
ties of the Ising model, the Baxter- Wu model and the Turban 
model with n = 3 and n = 4. x: bulk scaling dimension, 
z: dynamical scaling exponent, A: autocorrelation exponent, 
9: initial slip exponent. For n = 4 the phase transition is of 
first order: discontinuity fixed-point value, stretched 
exponential decay. Lower part: nonequilibrium critical ex- 
ponents of the two-dimensional Ising model starting with an 
initial state corresponding to the critical state of the Baxter- 
Wu model and the Turban model with n = 3 and n = 4. 
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FIG. 1: Typical initial state of the BW-model (left) and the 
n = 3 Turban model (right) with L — 120 having a magne- 
tization rrii — 0. The initial state is a mixture of the pure 
phases. 



perature as the initialization but the couplings between 
the spins are changed. For a given starting configuration 
the relaxation is repeated with typically a few hundred 
independent sets of random numbers. 



at the transition point, m c , through a stretched exponen- 
tial decay. On the other hand relaxation of the magne- 
tization starting with an uncorrelated initial state with 
a small magnetization, m„ has been shown to approach 
zero with an asymptotic power law time dependence, thus 
from a nonequilibrium point of view the transition is con- 
tinuous. 



III. RELAXATION FROM SECOND ORDER 
TRANSITION POINTS 
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In the following we present results of the relaxation 
of the critical Ising model where we start from a critical 
state of the BW model and of the n = 3 Turban model. 
The two models have the same static critical exponents 
and, interestingly, the decay of the critical correlations 
is the same as in the critical Ising model. We start to 
present first the results for the BW model. 



A. Relaxation from a BW critical state 

In the actual calculations we used finite systems com- 
posed of L x L spins, with L ranging from 60 to 240. 
We set the temperature to the critical value given in 
Eq. (j4j and then let the system evolve under the Hamilto- 
nian Hbw, see Eq.([3|), using the cluster-flip Monte-Carlo 
algorithm 36 . After equilibrium is reached the critical 
BW-states are extracted and sorted by magnetizations. 
For each magnetization, rrii, we have selected 1000 inde- 
pendent starting states. A typical starting configuration 
for rrii — is shown in the left panel of Fig. [TJ 

After having selected the initial state, we subject the 
system to the heat-bath dynamics of the critical Ising 
model, i.e. the relaxation is performed at the same tem- 



FIG. 2: (Color online) (a) Relaxation of the magnetization 
starting with a critical BW state with different initial magne- 
tizations rrii (from top to bottom, rrn — 0.33, 0.16, 0.08, 0.04, 
and 0). The data shown have been obtained for L = 60. (b) 
Relaxation of the magnetization starting with a critical BW 
state with initial magnetization rrii — in a log-log scale for 
different system sizes. 

Fig. [2] shows the time dependence of the magnetiza- 
tion for different starting magnetizations m^. It is seen 
in Fig. [2k that for small times the magnetization is 
monotonously decreasing, yielding for njj < 0.08 nega- 
tive values. The asymptotic dependence of m(t) is stud- 
ied in more details in Fig. ^jp for m t = where we 
plot the magnetization vs. time in a log-log scale. Ac- 
cording to this figure the absolute value of the magne- 
tization has a power-law dependence. For times shorter 
than some size-dependent time t\{L) there is an effec- 
tive exponent, G\ = 0.13(1), which turns to a different 
value, 9 — 0.18(1), for t > ti(L). The available time- 
scale is restricted by the finite size of the system, through 
t m r*> L z , z being the dynamical exponent of the critical 
Ising model, see Table HI The cross-over behavior seen in 
Fig. [2b is attributed to the structure of the initial state. 
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As seen in Fig. QJi the initial state can be described as 
a composition of the four different pure ordered phases 
of the BW model, one with magnetization, m = 1, and 
three others each with magnetization, m = —1/3. In 
the early time steps these pure phases relax and the re- 
laxation of the mixture of the m = —1/3 phases has an 
effective exponent, 9\ = 0.13(1). This assumption is in 
accordance with the effective exponent measured in Fig. 
[2^ for m,i — 1/3. Then for t > ti(L) the remains of the 
pure m = —1/3 phases are dissolved and we are in the 
true asymptotic regime. 

Next we consider the nonequilibrium autocorrelation 
function, G(t), calculated from a critical BW state with 
vanishing magnetization, m, = 0. The autocorrelation 
function is thereby defined by 



G(t) 



V- 



(7) 



where we average both over initial states and different 
realizations of the noise. According to the numerical re- 
sults shown in Fig. [3^, the autocorrelation function in 
a finite system of linear size L is well described by the 
functional form 



G(t) =At- x / z eM-t/r L ) , 



(8) 



where the characteristic time, tl, is a monotonously in- 
creasing function of L, so that in the thermodynamic 
limit G(t) has a power law dependence. We tried to fit 
the measured tl in the form t£ ~ and for small L 
the obtained exponent £ is compatible with the dynami- 
cal exponent of the Ising model, as listed in Table HI For 
larger L, however, when the characteristic time becomes 
larger than t\, the measured £ is found to decrease below 
1.5. A possible explanation of this behavior is that the 
relevant length scale in the problem, the typical size of 
the clusters of pure phases, £, see Fig. QJi is smaller than 
L. 

In order to deduce the autocorrelation exponent, A, 
from the numerical data we have calculated the cor- 
rected autocorrelation function, G(t) — G(t) exp(t/rr,), 
in which we used the estimate for tj, calculated from Fig. 
[3k.. G(t) is plotted in Fig. [5b in a double-logarithmic 
scale for different system sizes. Clearly, the curves are 
seen to approach an asymptotic curve for large L. 

Here, as for the relaxation process, one can observe 
an early time period in which case the effective expo- 
nent is Ai/z « 0.14(1). In the true asymptotic range the 
measured autocorrelation exponent is somewhat larger, 
X/z = 0.18(1). Interestingly, this value is much smaller 
than the value 0.74 obtained when starting from a fully 
disordered initial state. 



B. Relaxation from an n = 3 critical state 

For the case of the n — 3 Turban model we proceed as 
for the BW case. Using Glauber dynamics, we generate 
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FIG. 3: (Color online) (a) Autocorrelation function start- 
ing from a critical BW initial state for different system sizes. 
In the semi-log plot the asymptotic slope is proportional to 
the inverse characteristic time t(L). (b) The corrected au- 
tocorrelation function, G(t) = G(t) exp(i/ri), in a double- 
logarithmic scale for the same system sizes. The slope of the 
curves corresponds to the exponent X/z. 



equilibrium critical states of the Hamiltonian © with 
three-spin interactions. After equilibrium is reached the 
initial states are sorted by magnetization (see Fig. [TJd for 
an example). These initial states are then subjected to 
the heat-bath dynamics of the critical Ising model. The 
linear sizes of the system used here are the same as for 
the BW-relaxation described in Sec. IIII Al 

Starting with initial states with m, = the magne- 
tization relaxes to negative values. To obtain a quali- 
tative estimate of the initial slip exponent we plot in 
Fig. 2k the absolute magnetization vs. time in a log-log 
scale. From the asymptotic slope of this curve we obtain 
9 = 0.18(1), which is within error bars the value obtained 
from a BW initial state in Sec. IIII Al 



-l.o 

-1.5 
-2.0 



-3.0 



-4.0 



1 1 1 


11 Us- 4 














— jy * 
// / 






ff/ 
f* 






It 
— ft 






// 






-ft 
h 








-- L=60 






- L=120 






— L=240 




- (a) 
, 1 , 


1,1,1 





+ 

g 

"ic 



4 6 

In t 



10 




-2.0 - 



-2.5 



FIG. 4: (Color online) (a) Relaxation of the magnetization 
starting with a critical n = 3 state with initial magnetization 
rrii — in a log-log scale for different system sizes, (b) The 
corrected autocorrelation function, G(t) = G(t) exp(t/Ti), in 
a double-logarithmic scale for the n — 3 model. The slope of 
the curves corresponds to the exponent X/z. 

For the autocorrelation function the same ansatz as 
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given in Eq.© works in this case too, and the charac- 
teristic time, tl, has a similar size dependence as for the 
BW-model. Having estimated tl for each size we have 
calculated the corrected autocorrelation function, G(t), 
which is plotted in Fig. [4)d in a log-log scale. From the 
asymptotic slope of these curves we have estimated the 
autocorrelation exponent as X/z = 0.165(10), which is 
consistent , within the error of the calculation, with the 
value obtained from BW initial states in Sec. IIII Al 



IV. RELAXATION FROM A FIRST-ORDER 
TRANSITION POINT 

The Turban model with four-spin product interaction, 
n = 4, has an eight-fold degenerate state in the ordered 
phase, and this degeneracy is lifted at the first-order tran- 
sition point at T = T c . The transition point is char- 
acterized by the phase coexistence of the paramagnetic 
phase and the ordered phase, yielding a state that has a 
completely different structure as the equilibrium critical 
states of the n = 3 model or the BW model considered in 
the previous section. Since the order-parameter is non- 
zero even at T — T c , we consider in the following only the 
nonequilibrium autocorrelation function. Starting with 
an initial state of the n = 4 model at T c we have used 
heat-bath dynamics of the critical Ising model. 
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FIG. 5: (Color online) (a) Autocorrelation function starting 
from a first-order transition equilibrium state of the Turban 
model with four-spin product interaction for different system 
sizes. In the semi-log plot the asymptotic slope is proportional 
to the inverse characteristic time, t(L). (b) The corrected 
autocorrelation function, G(t) = G(t) exp(t/ri,), in a double- 
logarithmic scale for the same system sizes. The slope of the 
curves corresponds to the exponent X/z. 

The nonequilibrium autocorrelation function of this 
system is shown in Fig. [5^. The decay of the autocor- 
relation is again well described by the same functional 
form (jSJ that we have already observed for the two cases 
discussed before. From the system size dependent char- 
acteristic times tl, we infer a value z « 2.1 for the dy- 
namical exponent, in good agreement with the literature 
value 2.17 for the two-dimensional critical Ising model. 



The corrected autocorrelation G(t) = G(t) exp(i/TL) dis- 
plays a power-law decay as shown in Fig. [5Jd, but with an 
exponent X/z — 0.475(10) that is much larger then the 
value we have obtained when starting from the critical 
BW or n = 3 Turban model. Interestingly, the value of 
X/z is compatible with 1/z, yielding the value A = 1 for 
the autocorrelation exponent. 

In the following we present a simple scaling argument 
which could explain the origin of this numerical finding. 
First we note that the initial state at t = 0, being the 
state of a system at a first order transition point, is a 
mixture of clusters of the pure phases and the typical 
size of the clusters, £, is finite and given by the corre- 
lation length at the first order transition point. Now, 
for t > 0, we change the form of the interaction and 
let relax the system during which correlated domains of 
typical size £(t) ~ i 1 / 2 are created. In the correlated vol- 
ume, V(i) ~ £(t) d , due to random fluctuations, a given 
initial pure phase has an excess volume AV(t) ~ £(t) d / 2 
(compared to the average). Consequently the (corrected) 
autocorrelation decays as G(t) ~ AV(t)/V(t) ~ t~ d / 2z , 
so that A = d/2, in agreement with the numerical results. 



V. DISCUSSION 

In this paper we have considered a relaxation problem 
during which the form (and the symmetry) of the interac- 
tion between the particles is suddenly changed, whereas 
the temperature of the system is kept constant. In par- 
ticular we studied such cases when the system before and 
after the quench is at a critical temperature, but never- 
theless belongs to different universality classes. We have 
then studied well-known questions in nonequilibrium dy- 
namics, such as the relaxation of the magnetization and 
the decay of the autocorrelation function. It is known 37 
that nonequilibrium critical dynamics at time i, after a 
quench at t = from a state with Tj = oo, is analogous 
to the static critical behavior of a semi-infinite system at 
a distance y from a free surface located at y = 0. The 
analogous static critical problem to our dynamical prob- 
lem here is the interface critical behavior at a distance, 
y, from a straight interface, which separates two cou- 
pled semi-infinite critical systems which belong to differ- 
ent universality classe o 39 i 40 . According to our numerical 
calculations the nonequilibrium critical behavior in our 
problem is the result of the interplay and competition 
between the critical fluctuations of the two systems. 

The specific problem we studied is the Ising model on 
the square lattice, with different types of multispin inter- 
actions in the initial states, but with nearest neighbor in- 
teractions after the quench. Two initial models, the BW 
model and the n = 3 Turban model, belong to the same 
static universality class. Interestingly, the magnetization 
scaling dimension of these models, x = 1/8, coincides 
with that of the normal Ising model. Nevertheless due to 
the quench, the change of the symmetry of the interac- 
tion has a strong effect on the nonequilibrium dynamics. 
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Our main observation is that nonequilibrium dynamics 
has the same asymptotic behavior when we start from 
initial states of different models that are in the same 
(static) universality class. For the specific problem we 
considered here, the nonequilibrium exponents are found 
to satisfy the relation 8 — X/z within the error of the 
calculation. 

The third initial model we considered, the n = 4 Tur- 
ban model, has a first-order phase transition, thus the 
structure and the topology of the initial state is com- 
pletely different from the previously discussed cases. Also 
the asymptotic behavior of the nonequilibrium autocorre- 
lation function is described by a different exponent A = f . 
This exponent is expected to have the value A = d/2 and 
therefore to depend only on the dimension of the sys- 
tem, but to be universal otherwise, i.e. not dependent 



on the type of the initial state, provided it corresponds 
to a first-order transition point. 
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